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ABSTRACT 

The dynamical evolution of stellar clusters is driven to a large extent by their en- 
vironment. Several studies so far have considered the effect of tidal fields and their 
variations, such as, e.g., from giant molecular clouds, galactic discs, or spiral arms. In 
this paper we will concentrate on a tidal field whose effects on star clusters have not 
yet been studied, namely that of bars. We present a set of direct iV-body simulations 
of star clusters moving in an analytic potential representing a barred galaxy. We com- 
pare the evolution of the clusters moving both on different planar periodic orbits in 
the barred potential and on circular orbits in a potential obtained by axisymmetrising 
its mass distribution. We show that both the shape of the underlying orbit and its 
stability have strong impact on the cluster evolution as well as the morphology and 
orientation of the tidal tails and the sub-structures therein. We find that the disso- 
lution time-scale of the cluster in our simulations is mainly determined by the tidal 
forcing along the orbit and, for a given tidal forcing, only very little by the exact shape 
of the gravitational potential in which the cluster is moving. 

Key words: methods: iV-body simulations - galaxies: star clusters, - Galaxy: open 
clusters and associations: general 



1 INTRODUCTION 

While the internal dynamical evolution of star clusters is 
mainly driven by relaxation processes, the global galactic 
environment in which they are hosted can drive their evo- 
lution externally. Different aspects of star cluster evolution 
in the gravitational field of galaxies have been studied nu- 
merously using numerical simulations. Recent direct iV-body 
simulations include the effects of the Galactic tidal field 
(Kiip per et al. 2008; lErnst et al.1 12009ft and dynamica l fric- 
tion t erms for clusters close to the galactic centre (|just et al.l 
l201lh . From such kind of simulations the picture emerged 
that star clusters steadily loose stars through the Lagrangian 
points in the system and form long tidal tails. In axisym- 
metric potentials the stars in these tails are found to follow 
orbits which can be approximated by epicycles, resulting 
in a cycloidal motion of stars along the tails when seen in 
the ti dal tail coordinate system as introduced by I Ernst et al.l 
( 2009). This results in the formation of sub-structures within 
the tidal tails in form of clumps, often referred to as epicyclic 



* E-mail: iberent@ari.uni-heidelberg.de (IB) 
f E-mail: lia@oamp.fr (EA) 



over- densities. Numerical simulations have also shown that 
such structures persist also in tails of star clusters on non- 
circular orbits, however, their local density and po sition with 
respect to the cluster in t his case vary with time (|just et al.l 
l201ll : lKiipper et alJliuToh . 

A detailed understanding of the formation, morphology 
and dynamics of tidal tails resulting from tidally disrupted 
dwarf galaxies and star clusters in galaxies are of great in- 
terest to observations e s pecially in the Milky Way (e.g., 
ignkirchen et al.1 l200ll : iBelokurov et al.1 120061: IGrillmairl 
2009 , etc.), but also in other galaxie s (e.g. JPeng et al.ll2002l : 
Martmez-Delgado et aI1l20od I20T0L etc.). Tidal streams are 
of particular interest for tracing the gravitational potential 
by rec onstruction the orbits of the clusters, as e.g., for Palo- 
mar 5 lOdenkirchen et al. 200lh. The underlying assumption 
for this meth od is that the tidal t ails accurately trace the 
cluster orbit ([Varghese et alJlioilh . 

The galactic potentials used in numerical simulations 
such as the ones mentioned above are usually confined 
to sphericity or axisymmetry for simplicity. However, disc 
galaxies such as the Milky Way show a variety of morpho- 
logical and dynamical sub-structures that are known to af- 
fect the evolution of star clusters. E.g., the effect of inter- 
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actions with giant molecular clouds and spiral arm passages 
on the cluster dis soluti on time-scal e has been studied by 
iGieles et all (|2006h and lGieles et al.1 (|200/t ). respectively. 

In this work we focus on the evolution of star clus- 
ters in the time-dependent gravitational potential of barred 
disc galaxies. Stellar bars are known to be the main inter- 
nal driver of disc galaxy evolution and are strong non-linear 
perturbations to the gravitational potential of disc galaxies, 
thus having strong effects on the underlying orbital struc- 
ture in both 2-d and 3-d. We present a set of direct iV-body 
simulations of star clusters moving on planar periodic or- 
bits in barred potentials and on circular orbits in (near-) 
axisymmetric potentials. We particularly focus on the for- 
mation and evolution of tidal tails in both scenarios. 

The paper is structured as follows: In Section [2] we de- 
scribe the numerical methods and initial conditions used in 
our iV-body simulations. Sect ion [3] then describes the results. 
We give conclusions in Section [4] 



2 NUMERICAL METHODS AND INITIAL 
CONDITIONS 

2.1 Direct iV-body code using GPU clusters 

We use a modified versi on of the publicly available direct 
iV-body code (^-Grape ([Harfst et al.ll2007l ). It provides hi- 
erarchical particle time-st eps and a fourth o rder Hermite 
integration scheme (e.g., Mak ino et alJ Il992l ). Originally, 
this code is written for high-performance computing clus- 
ters equipped with the special-purpose h ardware Grape-6 
(|Makino et al.ll2003l : iFukushige et al.ll2005h to accelerate the 
gravitational force calculation. It has been used in the past 
for simula tions of dense stella r systems such as galactic nu- 
clei fe.g.. |Berczik et all 120051 . 120061 : iBerentzen et al.1 120091 : 
iKharchenko et al.l l2QQ9 ) . In our current version, the code has 
been adjusted to run on Gr aphic Processing Un its (Gpus) 
using the Sapporo library ([Gaburov et al.l 1200^ 1 which al- 
lows to use the Gpu cards like the special-purpose hardware 
Grape. 

There are alternative direct iV-body codes freely avail- 
able with more advanced features, such as regularisation 
methods for an accurate integration of binaries, stellar evo- 
lutionary tracks for single and binary stars, etc. However, 
the present work is the first to study star cluster evolution 
in barred galaxies and therefore we restrict ourselves to rel- 
atively simple models. Simulations with higher resolution, 
particle mass functions and increasingly more input physics 
will be studied in the future. 

In order to prevent the formation of stellar binaries or 
multiples in our simulations, we adopt a small gravitational 
(Plummer) softening in the force calculation. As a trade- 
off we study the evolution of the star clusters in the pre 
core-collapse phase only, where binaries and multiples are of 
less dynamical importance. As mentioned before, we limit 
ourselves in this work to iV-body models with equal mass 
particles. This prevents some p ossible dynamical e ffects such 
as mass segregation (see, e.g., iKhalisi et alll2007l ). We also 
neglect the dynamical friction between the star cluster and 
the galaxy, which generally becomes important close to the 
galactic centre. We discuss the possible effects of this limi- 
tation when necessary and leave simulations including this 
extra physics for the future. 



In the next subsection we describe the analytic galaxy 
model that we have implemented in the ip- Grape as an ex- 
ternal potential and force field. 

2.2 Analytic barred galaxy model 

We use an analytic three-component galaxy model which is 
composed of an axisymmetric disc, a spherically symmetric 
bulge/halo and a triaxial bar component. The disc poten- 
tial <I>d in our model is rep resented by a three dimensional 
Miyamoto- Nagai potential ([Miyamoto &; Nagalll975h of the 
form 



$mn(x,2/,z) 



GMq 



: + y 2 + (A + \[WT, 
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where G is the gravitational constant, Mo is the total mass 
of the component, and A and B are the radial and vertical 
scale-length, respectively. 

The spherically symmetric bulge/halo potential <I>b of 
our galaxy mod el is represented by a Plummer potential 
(|Plummerl 1 19 111 ), i.e., also described by Eq. [T] with A — 
0. For the bar potential ^ba r we use a th ree-dimensional 
Ferrers' ellipsoidal potential (|Ferrerdfl877h which is based 
on a density distribution of the form 
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and a b ar > ^bar > Cbar , (3) 



where Mbar is the mass of the bar component and abar, 
&bar and Cbar are its semi-principal axes. Note that in our 
description the bar major axis is aligned with the x-axis of 
our coordinate system. The full analytic expressions of the 
gravitati onal poten t ial <I>b ar and corresponding forces can be 
found in IPfennigerl ([ 19841 ) . 

The fourth order Hermite integration scheme used in 
our code requires both the particle acceleration a^ and its 
time derivative (also called 'jerk'). We integrate the evo- 
lution of the star cluster in a reference frame which is co- 
rotating with the bar, with a constant pattern speed Qbar- 
The resulting equations of motion for a particle are then 
given by, e.g., iBinnev &; Tremainel ([ 19841 ): 

tti = Od - V$ g al - 2 (ffcbar X f i) - f&bar X (f&bar X V l ) , (4) 

where a c \ is due to the self- gravity of the stellar cluster, 
^gai = + + ^bar and the last two terms are the 
Coriolis and centrifugal force, respectively. The jerk is then 
given as: 



: J el - ~Q t V$gal - 2 (ffcbar X tt t 



f&bar X (f&bar X V i) 



For the iV-body realisation of the star clusters we use 
distribution fu nctions following a Plumme r density p rofile 
(|Plummedll9lTh or a King density profile (|Kindll966h . We 
found no qualitative differences in the evolution between 



Table 1. Galaxy model parameter in model units. 



Component 


M [xlO 4 ] A [xlO 3 ] 


B [xlO 3 ] 


disc 


16.4 3.0 


1.0 


bulge/halo 


1.6 


0.4 



Mo is the mass of the components given in the first column. A and 
B are the radial and vertical scaling parameters of the Miyamoto- 
Nagai potential. 



Table 2. Bar model parameters in model units. 



Component 


M har 


a bar ^bar c bar 


bar 


2 x 10 4 


6 x 10 3 1.5 x 10 3 0.6 x 10 3 



Mbar is the mass of the bar component. (Zbar, ^bar and Cbar are 
the bar's semi-principal axes. 

the two types of models in our set of simulations. Since our 
simulations with King models have been run with a higher 
particle number and also have been carried out to longer 
simulation times we limit ourselves in this publication to 
those models. To set up the (non-rotating) isotropic King 
model for a given core radius r c and central concentration 
Wo we use Walter Dehnen's mk king routine from the open 
source Nemo software package (|Teubenlll995t ). 

All simulations presented in this work have been carried 
out on the dedicated Gpu cluster Kolob at the University 
of Heidelberg and at the Gpu facilities at the Laboratoire 
d'Astrophysique de Marseille. 

2.3 Parameter and units 

We set the gravitational constant G to unity and the units 
of mass and length to M u = 10 6 Mq and R u = 1 pc, re- 
spectively. With this choice, the resulting units of time and 
velocity become r u = 15 x 10 3 yr and v u — 65.6 kms -1 . Al- 
though the iV-body models are intrinsically scale-free, one 
should bear in mind that both the galaxy and the star clus- 
ter must be scaled simultaneously when scaling to different 
units, e.g., physical units. Especially, a rescaling of the mod- 
els will also change the bar pattern speed. Throughout this 
work we stay with the model units and refer to physical units 
only when necessary. 

Our choice of model parameter is base d on the we ll 
studied galaxy model described in ISkokos et al 
The corresponding mass distribution in this model results 
in a rotation curve as typically observed for disc galax- 
ies. The parameter converted to our set of units are given 
in Tables [1] and [2] for the galaxy and for the bar, respec- 
tively. The pattern speed of the bar is correspondingly set 
to Qbar = 7.64 x 10 -4 rad r" 1 with the bar rotating counter- 
clockwise in the (non-rotating) lab-frame. The rotation pe- 
riod of the bar corresponds to Tb ar ~ 8224 r u . 

The parameters of the star cluster (s) are summarised in 
Table [3 The half-mass radius of our (isolated) 50000 (50k) 
particle King model corresponds to rh ~ 4 and the tidal 
radius to about rt ~ 34. The total energy of the iV-body 
system is about E c \ — — 2 x 10 -3 . 
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Table 3. Star cluster parameters in model units. 



N 


Mtot 


W r c e 


50000 


0.2 


7.0 1.0 0.001 



N is the particle number used in the iV-body model and Mtot is 
the total cluster mass. Wo and r c are the concentration parameter 
and the core radius of the King profile, respectively, and e is the 
gravitational softening length used in the simulations. 
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Figure 1. Characteristic diagram of our barred potential. We 
plot the y-intercepts (yo) of the main planar prograde periodic 
orbits as a function of Jacobi energy Ej. The black line is the 
zero- velocity curve (ZVC). The main periodic 2-d families are 
shown in different colours as given in the legend. We mark the 
orbits on which we evolve the star clusters, indicating the stability 
by filled (stable) and open circles (unstable). 

2.4 Planar periodic orbits 

In this work we place the star cluster in the mid-plane of the 
galactic disc on an orbit which is periodic in the co-rotating 
reference frame of the bar. These orbits are the backbone of 
stellar bars (in 2-d) because they can confine regular regions 
(invariant tori) in phase-space around them. Orbits located 
in these regions show similar shapes and characteristics as 
their underlying parent periodic orbit. Thus we would expect 
a similar evolution of the star cluster when placed on more 
general orbits. 

In Fig. [T] we show the characteristic diagra m for the 
main planar orbit families (see, e.g., iBinnev &; Tremaind 
Il984h in the barred potential. We plot the y-axis intercepts 
(yo, with y < 0) of the periodic orbits as a function of their 
Jacobi integral, Ej, which is given as Ej = E — JQbar- Here, 
E and J are, respectively, the specific energy and the angular 
momentum of the orbit. The Jacobi integral is a conserved 
quantity in the frame corotating with the bafl 

The green line shows the so-called x\ 

orbit family, following t he nomenclature of 

IContopoulos h Papavannopoulosl (|l980h . These orbits 
are elongated along the bar major axis and close after one 
revolution in the bar potential and two radial oscillations 

1 Due to the self-gravity of the TV-body system the Jacobi integral 
is not a conserved quantity in our simulations. 
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Figure 2. 2-d planar periodic orbits in the barred and the axisymmetric potential. The corresponding orbit name, Jacobi integral Ej 
(for the barred model) and family is given in each frame. The colour coding corresponds to the strength of the tidal forcing along the 
orbit. See main text (Section I3.4|) for details. Iso-potential contours of the effective potential are plotted with grey lines. 
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(2/1). Close to the co-rotation radius the x\ turns over to 
the four-periodic (4/1) orbit branch (cyan line). The blue 
and magenta line in Fig. [T] represent the so-called X2 and 
X3 orbit families, respectively. These orbits are a class of 
2/1 orbits which are elongated perpendicular to the bar 
major axis and are a signature of the presence of the inner 
Lindblad resonance(s). 

The set of orbits on which we place the star cluster in 
our simulations are marked and labelled in Fig. [T] where 
we indicate the stability using different symbols. A detailed 
discussion of the orbit s and their stabilit y in the potential 
used here is given in ISkokos et all (|20Q2ah . 

For comparison, we also run a second set of simulations 
for each of the selected orbits, putting the star cluster on a 
circular orbit about the galactic centre using an (nearly) ax- 
isymmetric potential. For the latter, we maintain the mass 
of the bar, but set the parameter a in Eq. [3] to 2.7 x 10 3 , 
which is the mean value of the three axes as given in Ta- 
ble El We then set b = 0.99 a and c = 0.90 a due to the 
required inequality of the three parameter (see Eq. [3J. The 
axisymmetric galaxy components remain unchanged. 

The radius of the circular orbits is determined as fol- 
lows: we calculate the orbit- averaged radial acceleration or 
on a test-particle in the barred potential over one orbital 
period T or b for each of the periodic orbits marked in Fig[T] 
Using a simple bi-section root-finding algorithm, we then de- 
termine the radius of the corresponding circular orbit which 
has the same mean radial force. Information about both 
kinds of orbits (barred and circular) and about our set of 
simulations is given in Fig. [2] and Table [H 



3 RESULTS 

In this section we describe the evolution of the star 
clusters first in isolation and then when placed initially 
on a periodic orbit in our barred potential. We also 
compare the results with simulations of star clusters on 
circular orbits in the corresponding axisymmetric poten- 
tial. Many of the features described in this section are 
better seen in animation, than in plots. For this reason 
we include animations of all the models described here under 
http : //lam . oamp . f r/research/dynamique-des-galaxies/ 
scientif ic-results/star-cluster-evolution/. 



3.1 Simulation A - isolated star cluster 

Before studying the star cluster evolution in the gravita- 
tional field of our three-component galaxy model, we first 
run a simulation of the isolated cluster, i.e., without any 
external gravitational force. 

In Fig. [3] we show the Lagrange radii, i.e., spherical radii 
that contain a fixed fraction of mass, as a function of time for 
the isolated 50k King model (cf. Table [3]). We define the cen- 
tre of the cluster as the centre of mass of all particles within 
a sphere of r = 0.2 around the particle with the highest mass 



1 j 

<2r = jj- ^2 a j' e ~R,ji where slj is the test particles acceleration 



For a constant time-step At one gets Nj = T or ^/At + 1. 



Table 4. Simulations 



Orbit E 3 



family yo [lO 3 ] 



A 


Isolation 


- 


- 


- 


Bl 


-75.5 




0.57 


812.3 


B2 


-75.2 


circular 


0.46 


707.5 


CI 


-75.5 


XI 


0.10 


1041.4 




-Do. O 


circular 






Dl 


-54.0 


XI 


0.52 


3178.8 


D2 


-49.1 


circular 


1.57 


2111.0 


El 


-46.0 


XI 


0.99 


4654.5 


E2 


-40.6 


circular 


2.17 


2773.9 


Fl 


-43.1 


XI 


1.37 


6075.4 


F2 


-36.6 


circular 


2.55 


3262.1 


Gl 


-43.5 


XI 


2.59 


8008.6 


G2 


-32.2 


circular 


2.09 


3985.0 


HI 


-40.5 


4/1 


2.99 


14733.1 


H2 


-26.3 


circular 


3.99 


5167.0 


Jl 


-39.2 


4/1 


2.10 


14782.0 


J2 


-26.7 


circular 


3.92 


5083.5 



The first column gives the orbit name for each selected orbit. Ej 
is the Jacobi integral. The third column shows the corresponding 
orbit family, yo is the y-axis intercept of the orbit and T or b is its 
orbital period. 
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Figure 3. Lagrange radii as a function of time for the isolated 50k 
King model (simulation A). Red solid lines show mass fractions 
of 5, 10, 20, 30, . . . , 90 per cent of the total cluster mass and green 
solid lines are fractions of 91, 92, 93, . . . , 99 per cent (from bottom 
to top). The half-mass radius is plotted with a black solid line. 
The cyan solid line corresponds to the farthest gravitationally 
bound cluster particle, i.e., the particle with the largest distance 
to the cluster centre and a negative total energy. 



denisty. The latter has been c alculated using routi nes from 
the software package falcON (|Dehnenll2000[ l2002h . Due to 
close two- and few-body encounters between particles, its 
outer envelope slowly heats up and expands, while its inner 
regions simultaneously start to contract slowly. As can be 
seen from Fig. 02 the model does not undergo core collapse 
during the integration interval of 50000 time units, or some 
750 Myr - i.e., roughly 5 bar rotations. 

We find that all stellar particles remain gravitationally 
bound to the isolated cluster until the end of the simulation. 
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This is most likely due to the small gravitational softening 
which prevents particles to reach escape velocities during 
close stellar encounters. The introduction of a gravitational 
softening in our models slightly reduces the strength of close 
particle encounters and prevents the formation of bound bi- 
naries. Since our simulations end clearly before core-collapse 
the dynamical impact of binaries can safely be neglected. 

For this isolated King model the total energy E c \ of the 
cluster is a conserved quantity and can be used to check 
the accuracy of the numerical integration. We calculate the 
relative energy error and find that the energy is conserved 
to better than 5 x 10 ~~ 3 per cent over the full integration 
time. 



3.2 Periodic orbits - simulations B — J 

In this subsection we present a mainly qualitative descrip- 
tion of our simulations. We will use these results to quantify 
in the next subsection certain aspects of the cluster evolu- 
tion and the formation and evolution of the tidal structures. 



3.2.1 Simulations CI and C2 

The underlying orbit in simulation CI belongs to the x\ fam- 
ily (Figs. [T] and [2} and is dynamically stable. We place the 
cluster initially at pericentre at a distance from the galactic 
centre of about R — y^ — 100 (Table |4} in the mid-plane of 
the disc. The cluster is then launched from pericentre with 
an initial centre of mass velocity according to this x\ orbit. 
Since it is started relatively close to the centre of the galaxy, 
it is initially not in dynamical equilibrium - mainly due to 
the presence of strong tidal forces - and adjusts itself in the 
very early phases of evolution. This effect is visible in Fig. 2] 
in which we show the evolution of the Lagrange radii as a 
function of time for this simulation. 

After its first apocentre passage at t ~ 260 the cluster 
starts to form tidal tails which become clearly pronounced 
for the first time after about half an orbital period (see Ta- 
ble [4|). The leading and trailing tails are oriented along the 
underlying parent orbit and are initially confined to regions 
mainly inside and outside this orbit, separated by it. This 
is due to the fact that the particles are leaving the cluster - 
depending on their Jacobi integral - through the Lagrange 
points Li and L2, where L\ and L2 lie between the clus- 
ter and the gala ctic centre and on the opposite cluster side, 
respectively fsee lErnst et alj|2009l for an illustration). 

When the cluster approaches again the orbital apocen- 
tre, the stars forming the leading tidal tail start to slow down 
following their pericentre passage. This leads to a "compres- 
sion" first of the leading tidal tail, and consecutively of the 
cluster itself and of the trailing tidal tail as they approach 
the orbital apocentre and slow down as well. The compres- 
sion of the tidal tails is reflected in a decrease of the La- 
grange radii (Fig. which become minimal each time the 
cluster just passed the orbital apocentre. In this phase the 
trailing tail is still undergoing compression due to its decel- 
eration, while the leading tail starts to expand again. The 
expansion of the tidal tails is reflected in an increase of the 
Lagrange radii which reach their maxima when the cluster 
passes pericentre. 

The effect of compression and expansion of the tidal fea- 
tures close to apocentre and pericentre passage has also been 



10000 



1000 




0.01 



2000 4000 6000 8000 10000 
time [model units] 




0.01 



10000 



20000 30000 
time [model units] 



40000 



50000 



Figure 4. Lagrange radii of the cluster in simulation CI (cluster 
on an x\ orbit with Ej = —75.5) as a function of time. The 
grey vertical lines indicate the orbital period in the bar reference 
frame. The top and bottom panels show the time evolution up 
to 10000 and 50000 time units, respectively. The layout of both 
panels is the same as in Fig. [3] 



found in recent iV-body simulations of star cl usters moving 
on ec centric orbits in axisymmetric potentials ([Kiipper et al.l 
2010). The compression/expansion of tidal streams and the 
resulting asymmetry when close to apocentre should thus be 
inherent to clusters moving on non-circular orbits. The os- 
cillations found in the Lagrange radii appear with the same 
period as the one of the periodic orbit and are therefore 
linked to the variations of the tidal field along the trajec- 
tory of the cluster. 

As expected, the effect of periodic compression and ex- 
pansion is indeed not found in our corresponding simulation 
C2 in which the star cluster moves on a circular orbit in the 
axisymmetric potential. In this case the particles escape at 
an almost linear rate from the star cluster. For comparison 
with simulation CI we show in Fig. [5] the time evolution of 
the Lagrange radii for our model C2. The expansion of the 
different mass shells is relatively smooth and does not show 
any oscillations such as those found in model CI. This is 
due to the fact that the star cluster moves along the cir- 
cular orbit with almost constant azimuthal velocity and is 
exposed to a constant tidal field. Note that the radius of 
the farthest unbound cluster particle (see Fig. El cyan line) 
initially corresponds roughly to the initial tidal radius of the 
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Table 5. Periods. 



2000 4000 6000 8000 10000 

time [model units] 



on C1 with E 



75.5 



10000 



20000 30000 
time [model units] 



40000 



50000 



Figure 5. Lagrange radii for the simulation C2 (cluster on a 
circular orbit corresponding to orbit CI). Layout is the same as 
in Fig. [J] The sudden increase in the cyan line at around t = 18000 
results from a particle getting ejected from the cluster and thus 
is no more the last bound particle. 



cluster rt = 35, and slowly increases as the cluster is losing 
mass. 

The maximum value for the Lagrange radii is given in 
both simulations, CI and C2, by the separation of the two 
apocentres. We find in our run CI that the tidal tails fill 
about half the length of the x\ orbit after t — 4880, or 
roughly 5 orbital periods, and fills the full orbit at about 
t = 7800, or about 8 orbital periods. 

In both simulations, CI and C2, sub-structures form 
along the tidal tails in form of clumps (or epicyclic over- 
densit ies) as described and analysed i n previous publica- 
tions (jKiipper et al.l 120081 : lErnst et al.l 120091 : iKiipper et all 
2010). The position and strength of these clumps varies with 
the compression and expansion phases along the orbit. As 
the tidal tail reaches its maximum length at pericentre pas- 
sage of the cluster, the relative clump density increases, i.e., 
shows the maximal density contrast between clumps and the 
tidal tails. 

The periodic density variations of the tidal tails are vis- 
ible also in the oscillations of the outer Lagrange radii (e.g., 
those corresponding to mass fractions larger than 90 per 
cent), with a period estimated by eye to about half the or- 
bital period T or b (see Fig.SJ). To show that the density varia- 
tions are indeed correlated to the orbital period of the parent 
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-75.5 


x 2 


812.3 


411 ± 1 


2 


Dl 


-54.0 
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XI 
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» 2 


Gl 


-43.5 


Xl 
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» 2 


Fl 


-43.1 


Xl 


6075.4 


3034 ± 28 


» 2 


HI 


-40.5 


4 : 1 


14733.1 


3682 ± 41 


» 4 


Jl 


-39.2 


4 : 1 


14782.0 


3766 ± 43 


4 



The six columns, from first to last, give the orbit name, its Jacobi 
integral, the corresponding orbit family, its orbital period T or b, 
the period of the expansion/compression measured using the La- 
grange radius Lgg which contains 98 per cent of the initial cluster 
mass, respectively, and the ratio between the two periods. 



periodic orbit we use a Fast Fourier Transformation (FFT) 
to determine the oscillation period from the time series of 
the Lagrange radius L 98 , i.e., the radius which contains 98 
per cent of the initial cluster mass. For our run CI we find 
an oscillation period to be Tl = 529 =b 1, which in fact cor- 
responds to roughly half the orbital period 0.5T or b of the 
xi orbit. Our results for this analysis are summarised in Ta- 
ble [5] for all periodic orbits presented here. We find that at 
the end of the simulation (after 50000 time units) the cluster 
in run CI (barred potential) has lost about 79 per cent of 
its initial mass, while the same cluster in run C2 (circular 
orbit) has lost about 62 per cent of its initial mass. 



3.2.2 Simulations Bl and B2 

The parent periodic x\ orbit in the previous simulation CI is 
elongated along the major axis of the bar. We now describe 
the evolution of the cluster on an x 2 orbit with the same 
Jacobi integral of Ej = —75.5, but which - in contrast to x\ 
orbits - is oriented perpendicular to the bar major axis and 
is less eccentric. As before, this orbit is dynamically stable. 

Qualitatively, the evolution of the star cluster on the 
X2 orbit (run Bl) and on the corresponding circular orbit 
are found to be similar to those in simulations CI and C2, 
respectively. As in the previous simulations we find forma- 
tion of tidal tails and sub-structures therein in both kinds 
of models, i.e., in both the barred and the axisymmetric 
potentials. 

The effect of periodic compression/expansion of the 
tidal tails is less pronounced in run Bl (although clearly 
present) than in run CI since the underlying x 2 orbit is less 
eccentric. This is also reflected in the amplitude of the os- 
cillations (compare Fig. [6] left panel). Again, the successive 
compression and expansion of the tidal tails correlates with 
the orbital period of the x 2 orbit (see Table EJ). 

In simulation B2 (circular orbit) we note that the lead- 
ing tidal tail and the trailing tidal tail are confined to regions 
inside and outside the underlying circular orbit owing to 
the particles epicyclic motion. Due to the strong non-linear 
perturbation of the bar in our simulation Bl the (linear) 
epicyclic approximation does not apply to the motion of the 
particles in the barred potential. As result of this we find a 
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Figure 6. Lagrange radii for the simulation Bl (left panel; cluster on an X2 orbit with Ej = —75.5) and simulation B2 (right panel; 
corresponding circular orbit). The layout of both panels is the same as in Fig. [3] 



higher fraction of particles crossing the underlying periodic 
orbit as compared to the circular orbit B2. 

In simulation Bl the tidal tails fill one full orbit at about 
t ~ 15000. The mass loss of the initial cluster mass is about 
73 percent in run Bl and about 69 percent in run B2. In 
subsection 13.41 we will quantify and discuss this difference 
in terms of the tidal field for all our simulations. 



3.2.3 Simulations Dl and D2 

In contrast to the previously described simulations, orbit Dl 
shows pronounced loops close to its apocentr e (see Fig . 13, 
which is a clear signature for a strong bar (lAthanassoulal 
1992). The orbit is stable but lies in the close vicinity of 
a short instability strip in the x\ characteristic curve (see 

Fignj. 

We find that the formation and evolution of the tidal 
tails in simulations Dl and D2 is for several orbital peri- 
ods qualitatively similar to that of cases with lower Jacobi 
integrals. Initially the forming tidal tails follow closely the 
underlying periodic orbit. After several apocentre passages 
of the cluster, however, we find that the tidal tails become a 
bit more diffuse compared to the previously described sim- 
ulations. We attribute this effect to the dynamical stability 
of the periodic orbit. In the particular case of simulation Dl 
the underlying periodic orbit is dynamically stable. Parti- 
cles in the tidal tails which leave the cluster, however, enter 
regions in phase space which correspond to an unstable pe- 
riodic orbit and are thus not confined to the invariant tori 
of the initial parent orbit anymore. 

As before, we find oscillations in the Lagrange radii (not 
shown here) whose period correlates with the orbital period 
of the Dl orbit (see Table [5}. In Fig. [7] we show the cluster 
mass in run Dl as a function of time. Here we use two def- 
initions of the cluster members, namely, (a) particles which 
have never crossed a radius corresponding to twice the initial 
cut-off radius R C ut at any given time (Fig. [7] red line), and 
(b) all particles that reside within 2 x R cu t at any give time 
(Fig. [7] green line). The oscillations found in the latter case 
have approximately the same period as that derived from 
the Lagrange radius Lgg. A similar behaviour is also found 
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Figure 7. Mass loss from clusters in simulations Dl and D2. Red 
and blue line: Remaining cluster mass in Dl and D2, respectively, 
after removing all particles which have ever crossed the sphere 
around the cluster density centre with radius 2 x .Rcut- Green 
line: Mass of all particles in Dl that reside within this sphere at 
the current time. 



(not shown here) in the previously described runs Bl and 
CI. We note that the mass loss / m i(t) in all our simulations 
consists of two main contributions, i.e., an initial exponen- 
tial decay fi(t) continuing with a secular linear decay f2(t): 



/mi (t) = h{t) + f a {t) = (mo e" rai ( ) + (m 2 - m 3 t) . (6) 

We discuss this dependency in more detail in Sec. I3.4I At 
the end of simulation Dl about 52 percent of the initial 
cluster mass is gravitational unbound (see Fig.0), while for 
the corresponding circular orbit in simulation D2 about 23 
per cent of the original cluster mass is unbound. 



3.2.4 Simulations El and E2 

To test our hypothesis of the influence of the stability of 
the periodic orbits on washing out the tidal tails we chose 
now an unstable x\ orbit roughly centred on the instability 
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Figure 8. Lagrange radii for the simulation El (left panel; cluster on an x\ orbit with Ej = —46.0) and simulation E2 (right panel, 
cluster on the corresponding circular orbit). The layout of both panels is the same as in Fig. [3] 
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Figure 9. Star cluster in simulation Dl at time t = 24000. The 
particles are coloured according to the local mass density of the 
cluster. The underlying periodic is shown with a red line. The 
arrow points to the cluster position. Leading and trailing arms 
are labelled accordingly. 

strip in the characteristic curve between Ej « —48.4 and 
—43.6. As in the previous case, this x\ orbit shows loops 
close to its apocentres (see Fig. [2}. In Fig. [8] we show the 
evolution of the Lagrange radii for simulations El and E2. 
The initial evolution of the cluster and its tidal tails follow 
the underlying parent orbit. After a few apocentre passages, 
however, the structure of the trailing tidal tail becomes more 
diffusive and even shows multiple tail- like structures (e.g., 
compare Figs. [9] and I10p . Both the leading and the trailing 
tidal tails wash out after a few apocentre passages, while the 
cluster on the corresponding run E2 (circular orbit) shows 
nicely separated tidal tails. 

The clusters in El and E2 have lost by the end of the 
simulation about 30 and 17 per cent of the initial cluster 
mass, respectively. This supports the trend that a star clus- 
ter in barred potential loses more mass than the same cluster 
on the corresponding circular orbit in the axisymmetric case, 
when being evolved for the same time period. 

3.2.5 Simulations Fl and F2 

The x\ orbit in our simulation Fl has peak- like edges at 
apocentre. The orbit lies right at the knee (kink) of the x\ 
characteristic curve in a stable region, but just next to an 
unstable region (towards larger pericentres) . 



imulation E1 E,= -46.0 

t = 07810 




Figure 10. Star cluster in simulation El at time t = 7810. The 
particles are coloured according to the local mass density of the 
cluster. The underlying periodic orbit is shown with a red line. 



The formation of the tidal tails and their evolution is 
similar to that of the previous models. At t ~ 7300 the 
leading tidal tail is highly compressed into one leading clump 
with distance of AR ^150. The star cluster in run Fl loses 
about 11 per cent and in run F2 about 14 per cent of its 
initial mass. 



3.2.6 Simulations Gl and G2 

This orbit is a stable x\ orbit. The formation of the tidal 
tails is similar to the case of circular orbits. The leading tail 
slows down when approaching apocentre and the clumpy 
sub-structure is stronger in this phase. The location of the 
clumps within the tails changes with time. Indeed, as the 
clumps in the leading tail slow down when approaching 
apocentre their relative distance with respect to the clus- 
ter centre becomes shorter. At the end of the simulation the 
tidal tails extend roughly over the underlying periodic orbit. 
As in the previous runs the centre of mass of the cluster does 
not stay on the underlying orbit. The star cluster in simula- 
tions Gl and G2 loses about 19 and 9 per cent, respectively, 
of the initial cluster mass. 
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Figure 11. Axial ratios c/a (red line) and b/a (green line) derived from moments of inertia within a spherical shell with r = 5. . .25, 
as well as the corresponding triaxiality parameter # (blue line) as given in Eq. [7] We also show the angle between the clusters centre of 
mass velocity vector and the major axis unit vector based on the moments of inertia measured (a) in shells with r = 5 . . . 25 magenta 
line) and (b) with r = 100 — 150 (£ ; cyan line). Left panel: results from run Gl. Right panel: results from run G2. 



3.2.7 Simulations HI and H2 

The orbit in simulation HI is a 4/1 orbit, i.e., a four-periodic 
orbit closing after one rotation and four radial oscillations. 
It is an unstable orbit and has sharp pointy edges at its 
apocentres. As before, we witness the formation of tidal tails, 
which after apocentre passages do not follow the underlying 
orbit anymore. The evolution of the star cluster in this run 
further supports our hypothesis that the orbit stability in 
the global galactic potential determines the orientation and 
shape of the tidal tails and also how diffusive they are. The 
cluster in run H2, i.e., on the corresponding circular orbit, 
does not show any peculiarities as compared to the earlier 
simulations using circular orbits. The mass loss at the end 
of the simulation is about 8 per cent in run HI and about 4 
per cent in run H2. 



3.2.8 Simulations Jl and J2 

This orbit in simulation Jl is again an unstable 4/1 orbit. 
In contrast to the previous orbit HI, however, it has loops 
at its apocentres. The tidal tails which form in this model 
are initially aligned with the underlying periodic orbit. This 
changes after the first apocentre passage, when reaching the 
outer edge of the loop there. The tail rapidly swings around 
when the cluster runs through the loop and aligns again with 
the orbit after the cluster leaves the loop. After the third 
apocentre passage the tidal tails are less sharply denned and 
less aligned with the periodic orbit. As in the other cases 
with unstable orbits, the tidal tails slowly dissolve into a very 
diffusive morphology. We furthermore note that the cluster 
itself does not remain centred on the initial periodic orbit, 
but its density centre reaches separations up to AR ^100 
from the initial underlying periodic orbit. 

Both the mass loss rate of the cluster and the density 
in the tidal tails is relatively low. At the end of the run 
the tidal tails have filled about half the underlying orbit. In 
simulation Jl the star cluster loses about 11 per cent of its 
initial mass and about 5 per cent in run J2. 



3.3 Shape and orientation of the cluster 

In order to quantify the shape and orientation of both the 
cluster and the tidal tails formed in the simulations, we cal- 
culate the moments of inertia tensor of the cluster inside 
a spherical shell within 5 ^ r ^ 25 and determine the 
corresponding set of eigenvalues I a ,Ib,Ic, with I a ^ h ^ 
I c . For the analys is we use the EiSPACK software package 
(jSmith et al.l ll976). From the eigenvalues we then determine 
the axial ratios as b/a = \J (I c + I a — h) / {Ic + h — la) and 
c/a — \/(h + I a — Ic) I {Ic + h — I a) as well as the triaxial- 
ity parameter as denned, e.g., in Fran x et al. 

I (|l99lh . as 

tf=(l-b 2 /« 2 )/(l-c 2 /« 2 ) , (7) 

with a > b > c. In terms of the eigenvalues, Eq. [7] can also 
be rewritten as $ = (h — la)/ (I c — la)- This definition gives 
$ = 1 for a prolate object (a > b — c; cigar-like shape) 
and $ = for an oblate object (a = b > c; disc- like shape). 
An example of the time evolution of $ is shown in Fig. [11] 
(left panel; blue line). We find that star clusters on orbits in 
barred potentials and high energies Ej typically have smaller 
triaxiality parameter while clusters on orbits with smaller 
energy Ej have values of $ « 1. 

Furthermore, we determine the angles Q between the 
current centre of mass velocity vector of the cluster and 
the major axis inertia eigenvector determined from the mass 
distribution of the cluster within a spherical shell between r 
= 5 and 25. Similarly we define £ using the corresponding 
eigenvector for the tidal tails (measured between r = 100 and 
150). Fig. 1111 fright panel) shows the results for the circular 
orbit G2. Once the tidal tails have formed, the axial ratios 
of the cluster remain roughly constant. From the moments 
of inertia we find that the cluster is elongated perpendicular 
(i.e., £i ~ 0.57r) to its trajectory (due to the radial tidal 
force). The tidal tails follow the underlying orbit. Within 
the region r=100 and 150 (which is relatively close to the 
cluster compared to the total length of the tails) we get 
Co ~ 0.257T for our run G2. 

For orbits in barred potential (see Fig. [11] left panel) 
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Figure 12. Results from our simulations CI (top panel) and El 
(bottom panel). The layout is the same as in Fig. 1111 
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the semi-major axis of the cluster lies in the orbital plane but 
periodically changes its orientation within the orbital plane 
with respect to the cluster trajectory. The same is true for 
the tidal tails. The oscillations of ( Q have a larger amplitude 
than those of The evolution of the orientation of both 
the cluster and tidal tails is more regular as compared to 
orbits in the barred potential. We note that the amplitude 
of the oscillations in £ G are higher for orbits which have large 
variations in their radius of curvature (compare Fig. I12[) . 

In Fig. [13] we plot the time-averaged values of the axial 
ratios < c/a > and < b/a >, as well as of the triaxiality pa- 
rameter < i9 > , all as a function of the Jacobi integral Ej of 
the underlying periodic orbit. We find that star clusters at 
larger Ej roughly maintain their spherical symmetry, much 
more so than clusters on orbits with small Ej (i.e. clusters 
near the centre), which become considerably triaxial. Clus- 
ters nearer to the centre also have larger dispersions of their 
axial ratios. The flattening jc/a£ of the cluster follows the 
same trend for simulations in barred and in axisymmetric 
potentials. This is not true for the elongation jb/a£ in the 
orbital plane, which has a stronger dependence on Ej in 
barred than in non-barred cases. These effects are presum- 
ably the reflections of the variation of the tidal force along 
the corresponding periodic orbit. 



Figure 13. Time-averaged axial ratios < c/a > (top panel), 
< b/a > (middle panel) and triaxiality parameter < $ >. The 
filled and open symbols represent simulations with barred orbits 
and circular orbits, respectively. The corresponding small symbols 
show the standard deviation, of these quantities. 



3.4 Tidal forces 

To quantify the strength of the tidal forcing along the peri- 
odic orbits we calculate the eigenvalues (<?i, (72, #3) and the 
corresponding eigenvectors of the force tensoiu, as given by 



(8) 



dx i dx j 

We find that one of the eigenvectors is always parallel to 
the z-axis, i.e., oriented perpendicular to the orbital plane. 
In the following we call the corresponding eigenvalue q z and 
the two remaining eigenvalues q a and The eigenvectors 



3 This definition is s imila r to that of the tidal tens or used in, e.g., 
iRenaud et all (120081 ) and lPasetto fe Chiosil <2009h . 
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Figure 14. Cluster expansion rate based on a linear fit to the 
(logarithmic) Lagrange radii of the outer mass shells. Filled and 
open symbols indicate simulations in barred and unbarred poten- 
tials, respectively. 



Figure 16. Final mass loss (in per cent) in unbound material. 
Filled and open symbols indicate simulations in barred and un- 
barred potentials, respectively. 
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Figure 15. Mass loss parameters as a function of the average 
tidal forcing. Filled and open symbols represent simulations in 
barred and axisymmetric potentials, respectively. The parameter 
of the linear fit /2(£) = m2 — rri3 t are plotted with circles (7712) 
and triangles (7713). The parameter of the exponential fit /i(t) = 
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-mi t 



are plotted with diamonds (mo) and squares (mi). 



of q a and o> lie within the orbital plane. Positive eigenvalues 
are found for expansive flows, while negative ones represent 
for compressive flows. We measure the tidal forcing on the 
cluster by the quantity y/ q% + q% . In Fig [2] we use this quan- 
tity for the colour coding of the orbits. Note that the tidal 
forcing in barred potentials is strongest along the bars inter- 
mediate axis 6bar for our simulations with Ej ^ —43.1 and 
along the bars major axis aw otherwise. 

As can be seen in the figures showing Lagrange radii, the 
expansion rate of the different mass shells typically seems to 
be higher for orbits in the barred potential as compared to 
the circular ones. In Fig. [14] we show the rate at which the 
cluster is expanding/dissolving. We determine this rate by 
using a linear fit to the steeply rising part of the Lagrange 
radii of the outer mass shells. We then calculate the aver- 
age expansion rate < dr^/dt > for each run from about 5 
determined slopes. The top panel in Fig. [14] shows the ex- 
pansion rate as a function of the Jacobi integral Ej . We find 
a strong correlation between the expansion rate and the Ja- 
cobi integral. Equivalently, we also plot the expansion rate 
as a function of the orbit- aver agecQ tidal forcing (see bottom 
panel in Fig. I14|) . Generally, the expansion rate is found to 
be higher for clusters in barred potentials (at least for small 
Ej). 



4 See Sec. 13. II for details. 
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Figure 17. Distribution of Jacobi energy Ej for simulation Dl. The times from top to bottom correspond to 16, 16 ^4, 36, and 44 x /4 
orbital periods, respectively. The left panels show the projected particle distribution colour coded by Ej. The underlying periodic orbit 
is shown by the blue line and the position of the cluster is marked with a black dot. The right panels show the binned Ej distribution 
using the same colour coding as before. The Jacobi integral of the underlying periodic orbit is indicated by the vertical thick line. 



In all our simulations the clusters lose their mass (i.e., 
stars become gravitationally unbound) roughly linearly. For 
clusters starting close to the centre of the galaxy, we ob- 
serve an additional, exponentially decaying mass loss, as a 
consequence of the non-equilibrium initial conditions of the 
cluster in the strong tidal field of the galaxy. To quantify the 
mass loss of the clusters in our simulations, we fit / m i from 
our Eq. [6] to the mass as a function of time curves, such as 
shown in Fig. The results of these fits is shown in Fig. 1151 
For simulations with a weak tidal forcing we get a ^/-intercept 



of the linear contribution of rri2 = 0.2, which corresponds to 
the initial cluster mass. This shows that, in this case, there 
is no significant contribution from the exponential mass loss. 
For higher values of the tidal forcing (and the initial tidal 
perturbation), the exponential part fi(t) becomes more and 
more important. At the same time, the value of rri2 drops 
to roughly 0.1. The linear mass loss rate (7713), shown as tri- 
angles in Fig. 1151 increases with increasing tidal forces. The 
solid triangles (barred models) lie slightly higher than the 
open ones (axisymmetric potential). 



14 /. Berentzen and E. Athanassoula 



The exponential mass loss rate mi (squares) is roughly 
independent of the tidal forcing, but the amplitude mo (dia- 
monds) increases considerably with increasing tidal forcing. 
To summarise, it seems that if we compare clusters in barred 
and axisymmetric potentials that there is on average very 
little difference in the overall mass loss rate, if the average 
tidal forcing along the orbits is the same. 

In Fig. [16] we show the total mass loss of the cluster 
determined at the end of the simulations. We again find 
strong correlations between this mass loss and the Jacobi 
integrals of the orbits or the average tidal forcing. 

3.5 Energy considerations 

In Figures [T3 we show the evolution of the energy distribu- 
tion within the cluster for simulations Dl (barred potential). 
This is of particular interest concerning the interpretation 
of observational dat a on the virial state of a star cluster 
(|Kiipper et alJ fcoil). e.g., such as Palomar 13. 

Fig. [TT] shows the particle distribution in simulation 
Dl (left panels) and the corresponding energy distribution 
(right panels) , both with a colour scheme based on their Ja- 
cobi integral. The snapshots are chosen at times when the 
cluster is at apo- (second and fourth rows) and peri-centre 
(first and third rows). With the formation of the tidal tails 
and the accompanied dissolution of the main cluster body 
the energy distribution slowly forms a double peak structure 
roughly symmetric with respect to the Jacobi integral of the 
periodic orbit. Thus the whole distribution has a three peak 
structure: the two peaks corresponding to the tidal tails, plus 
the central peak of the initial cluster distribution. The latter 
can at times be relatively low, so that the central structure 
may look like a plateau, particularly with the logarithmic 
scale we are using. At pericentres the double peaks reach 
higher values than at apocentres. However, the relative con- 
trast between the peaks becomes much strongest at apoc- 
entres, so that they stand out clearest. This is due to the 
compression effect described in Sect. 13.2.11 Note that the 
central peak is not pronounced at pericentres, but stands 
out clearly at apocentres, where it is an order of magni- 
tude higher than at pericentre. Comparing these structures 
always at the same phase of the orbit (i.e. always at apocen- 
tre, or always at pericentre) we find that, as time increases, 
the height of the two side peaks increases, that of the central 
peak decreases and the contrast between the peaks becomes 
much sharper. This is due to the fact that more and more 
material leaves the cluster and moves to the tails, as seen in 
Sect. [3 

Note that the double peak structure mentioned above 
is also visible in simulations with circular orbits, but is nev- 
ertheless much less pronounced - due to the absence of the 
compression/expansion along the circular orbit. Further dis- 
cussion of the tails and of the corresponding Ej distributions 
will be given elsewhere (Athanassoula, Romero-Gomez and 
Berentzen, in prep.). 



4 SUMMARY AND CONCLUSIONS 

Stellar streams are fossils of the formation history of galax- 
ies and can be used to reconstruct, e.g., the shape of the 
gravitational potential of dark matter halos. Several issues 



such as the life-time of the clusters and their tidal streams 
however are not yet fully understood. The understanding 
of the morphology and the dynamics of the tails and their 
sub-structures is still in its infancy. Numerical A^-body simu- 
lations are a main key to address these questions. Many sim- 
ulations in the literature focus on the evolution of dense stel- 
lar systems in an external axisymmetric potential on both 
circular and eccentric orbits, respectively. For orbits close 
to the galactic plane, however, non-axisymmetric perturba- 
tions in the disk, such as stellar bars and/or spiral arms, are 
expected to significantly influence the dynamical evolution 
of star clusters. 

In this work we have presented a set of direct A"-body 
simulations of star clusters orbiting in the plane of a barred 
galaxy potential. The clusters are placed in the mid-plane 
of the galactic disc and launched on periodic orbits which 
have been selected from the main planar 2-d families. The 
results are then compared to the evolution of star clusters 
on circular orbits in an axisymmetric potential. 

We find that the star clusters dissolve due to the tidal 
force field of the galaxy - in both the axisymmetric as well 
as in the barred potential. While the tidal tails in the ax- 
isymmetric potential continuously grow in length, we find a 
periodic compression and expansion of the tails for star clus- 
ters in the barred potential. Both the length and the density 
of the tidal tails vary with time and are directly correlated 
to orbital period of the underlying parental orbit. We find 
sub-structures within the tidal tails in form of clumps (or 
epicyclic over-densities). The qualitative evol ution of these 
struct ures is in agreement with the results of IKiipper et al.l 
( 2010), who studied the evolution of star clusters on eccen- 
tric orbits in a static axisymmetric potential. 

We also studied the distribution of the Jacobi energy 
Ej of the cluster stars and found that it depends strongly 
on the location of the cluster along the orbit. It remains, 
nevertheless, roughly symmetrical with respect to the Jacobi 
energy of the cluster orbit. Two peaks develop, one on each 
side of this energy, and their relative amplitude depends on 
the location of the cluster along the orbit. At apocentres 
the distribution is tri-modal and the three peaks are very 
clearly separated. At pericentres, the distribution is much 
more extended and the central peak is not clearly delineated. 

We find that the mass loss of the cluster is mainly deter- 
mined by the orbit- averaged tidal forcing. In other words, 
the shape of the gravitational potential - at least in the 
case studied here - affects the dissolution time-scale only 
very little. On the other hand, the stellar bar has a strong 
effect on the shape and morphological evolution of the tidal 
tails: the stability of the orbit has a strong impact on how 
long the tidal tails and their sub-structures actually survive. 
This makes it more subtle to reconstruct the cluster orbits 
in observations and thus to probe that galactic potential in 
the disc plane. This could also be the case small galactic alti- 
tudes in which the 3- d shape of orbit s is strongly affected by 
the stellar bar (e.g. , [Pfenniger 1984; Be rentzen et al.l fl998; 
ISkokos et alJl2002aL lb nPatsis et al.ll2002h . 
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